BELTRAMI FORMS, AFFINE SURFACES AND THE 
SCHWARZ-CHRISTOFFEL FORMULA: A WORKED OUT 
EXAMPLE OF STRAIGHTENING 



ARNAUD CHERITAT 

Abstract. Consider the straightening </> of a Beltrami form that is constant 
on a square, with the corresponding ellipses having a vertical major axis, and 
null outside. A generalized Schwarz-Christoffel formula is used to express the 
inverse of (f). The formula is found by introducing an affine Riemann surface. 
This formula is used to draw on a computer the image of the square by (/!>, and 
practical aspects are discussed. The resulting shapes are shown for different 
values of the constant dilatation ratio of the ellipses (=major axis/minor axis). 
The limit when this ratio tends to infinity is surprising. A model of this limit 
is proposed, produced by an affine surface uniformization. 



In all this article, unless otherwise stated, the word afhne means complex-affine. 
The term "affine map" will alway refer to non constant afhne maps. The word 
surface will not refer to complex manifolds of complex dimension 2. 

1. Structure of the article 

Section [2] plays the role of an introduction. There, we present the problem that 
motivated this work. It involves a Beltrami form that is constant in a square and 
zero outside. What does the image of the square look like? What is the limit of this 
shape when the dilatation ratio tends to infinity? The author shows a few early 
computer experiments. 

In section [3] we explain how to express the inverse of the straightening of the 
Beltrami form using a generalized Schwarz-Christoffel formula. The problem of 
straightening the Beltrami form is transformed into the problem of conformally 
uniformizing to a punctured Riemann sphere some affine surface. We introduce 
a differential invariant characterizing the induced affine surface structure on the 
punctured Riemann sphere; it turns out to be a simple rational map. From this we 
deduce the formula. 

In section |4] we discuss practical aspects of numerically using the formula. 

In section [5] we study the limit of the formula as the dilatation ratio tends to 
infinity, and give a candidate affine surface whose uniformization would correspond 
to that limit. 

2. A TEASER 

In this section, we will present two questions, that originally motivated the work 
of the author, and show a few computer experiments. The tools to answer the sec- 
ond question are introduced in section [31 This will allow us to run better computer 
experiments in section [4.3.41 helping us to guess the answer. The same tools allow 
to give a precise candidate for the answer, and this is discussed in section [5l 

In subsection 1 2 . II we state the questions. In subsection 12. 21 we begin exploring it 
with computer experiments. In subsection 12.31 we give a few nai've guesses of the 
possible answer. 
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2.1. Straightening a square. Consider the square Sq in C defined by x + iy G Sq 
<J=4> X e [—1, 1] and y £ [—1, 1]. Consider the Beltrami form fJ-iz)^ with = 
outside Sq and fi{z) = i in Sq. The corresponding ellipses are circular outside Sq 
and vertical with ratio K = 2 in Sq. The straightenings of this Beltrami form are 
quasiconformal over C, and conformal outside the square. Let (j) : C —> C he the 
unique straightening satisfying the following normalization: 

(l){z) = z + + o(l). 
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Figure 1: The Beltrami form fj,, visualized as a field of ellipses. 

First question: what does the image of the square look like? 

One may give a rough guess, with the following heuristics: If one zooms on 
one corner, say the upper right one, /i tends to some Beltrami form i/{z)^ with 
fiz) = 1/3 on the lower left quadrant, and anywhere else. This Beltrami form is 
explicitly solvable, by the following method, illustrated on figure^ slit the complex 
plane along the vertical half-line L — {~oo,0]i. Define a map / :C\L^C\Lby 
f{x + iy) = X + in the lower left quadrant, and f{z) — z elsewhere. This map 
solves the Beltrami equation, but has a discontinuity along L: a point that crosses 
this boundary from the lower left quadrant to the lower right quadrant sees its image 
by / jump from iy/2 £ L to iy. The idea is then to define a Riemann surface iS 
obtained by gluing the lower left and the lower right quadrants along their common 
boundary L \ {0} by identifying z/2 to z. It is known that this Riemann surface is 
uniformized to C* by a branch of the map z i— s- z" for some well chosen complex a. 
Indeed in the coordinates w = logz, the problem is equivalent to gluing the band 
"Imz e [— 7r/2, 37r/2]" hy x — i7r/2 ^ x + 2>i'iT/2 — log 2, i.e. by the translation of 
vector 2i7r — log2. The job is done by for instance hy w ^ exp(2i7rw/(2i7r — log2)). 
So a = 2i7r/(2i7r - log 2) = 1/(1 - log(2)/2i7r). The image of the two half-lines 
composing the boundary of the lower left quadrant (and of any half-line through 
0) are logarithmic spirals. They turn quite slowly: to make one turn, the distance 
to needs to be divided by e'-^'^^ /'"s^ « 6 x 10^''. So we may expect the image of 
the square to look roughly like Figure [31 

Second question: what if one lets K tend to infinity? 

Since is conformal outside Sq and 4>{z) '--^ z at infinity, the capacity of 0(Sq) 
is equal to that of Sq, and univalent function theory implies that the image (^(Sq) 
stays contained in some disk B{0,R) independent of K. What else can be said? 

There is a not-so-well-known similar case with a very simple solution: replace 
Sq with the unit disk D, let /i be outside D and constant /i = a e (0, 1) within D, 
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Figure 2: Straightening the Beltrami form corresponding to horizontal ellipses with 
ratio 2 on the gray quadrant and circles outside. 




Figure 3: A sketch of what the image of the square is expected to look like, with 
the spirals exaggerated. It was hand drawn with a Vector Drawing Program. 

corresponding to vertical ellipses with ratio K = jz^- Then the straightening has 
the following form: if \z\ < 1 then (j){z) = z + az, if \z\ > 1 then (f){z) = z + a/z. The 
image ^(ID)) of the unit disk is a horizontal ellipse, the lengthes of its semi major 
and semi minor axes being I + a and 1 — o. When K — > +oo, i.e. when a — > 1, 
this ellipse tends to the segment [—2, 2] and converges uniformly. Its limit is z + z 
in D, which flattens D into [—2,2], and z + 1/z on C \ D, which is a conformal 
mapping from the exterior of the unit circle to the complement of [—2, 2]. 




Figure 4: Straightening a Beltrami form which is constant in the unit disk and zero 
outside. 

But for the square, it cannot be so sirnple. For one thing, its right side segment 
accounts for 1/4 of the harmonic measurqj of dSq with respect to oo. Therefore, its 

^Harmonic Measure has several meanings in mathematics, so let us precise the one we mean. 
Given a connected compact subset K with more than one point and with connected complement, 
there is a unique p > and a unique conformal bijection (f> from C \ -fC to C \ pB such that 
4>{z) — z tends to at infinity. The potential associated to K is the function G{z) = log |(/>(2;)|. 
It turns out that G{z) = log \z — u\dm{u) for a unique measure m with mass 1 and support 
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image by (j) must account for 1 /4 of the harmonic measure of 4>{dSq) with respect 
to oo. Thus it cannot tend to a point. Otherwise the capacity would have to tend 
to +00. 

2.2. First computer experiments. To get a better idea of what is going on, 

the author, which will be referred to as T in the rest of this section, decided to 
make computer experiments, and took a Partial Differential Equation approach. 
Indeed, we can consider (p as & solution of = 0, for a modified Laplacian A. 
The reader will sec from the formulae below that the modified version still satisfies 
Re (A0) = A(Re0) and Im(A(/)) = A(Im(/)), so solving At/) = for the complex 
valued function ^ amounts to independently solve the equation for its real and 
imaginary parts. 

Far enough from 0, (l){z) is close to z. So I took a big enough square, of side 
R, centered on 0, and tried to numerically solve the following Dirichlet problem on 
a fine enough grid dividing this big square: A(t> = with data = z on the 
boundary of the big square for some modified discrete Laplacian A. 

Usually, on a grid indexed by G 1? , given a function {i,f) 1— > m^j, the 

standard discrete Laplacian has expression Au : (i, j) 1— > Ui^\_j + Ui+x^j + + 
—4uij. The modified discrete Laplacian is defined as follows. Let z = Zij be 
the point on the grid corresponding to index {i,j). This point Zij can be outside Sq, 
inside Sq or on one of its four sides U, D, R, L, or on a corner. Let us introduce the 
following local coordinates that locally solve the Beltrami equation near the point 
{i,j): first we translate by —Zi^j to put z on the origin. Then we compose with a 
function that depends on where is z: 

• we compose with the identity if z is outside Sq; 

• with (x + iy) 1— > Kx + iy if z is inside Sq; 

• if 2; is on the left side of the square Sq, we compose with id on the left 
half-plane and Kx + iy on the right half-plane; 

• if 2; is on the right side, we take them the other way round; 

• if 2; is on the upper side, we compose with id in the upper half plane and 
X + iy/K in the lower half plane; 

• if 2; is on the lower side, we take them the other way round; 

• if 2; is on a corner. . . let us temporarily put corners aside. 

In each case, let us call I, r, u, d the length of the image, by this composition, of 
the segments from Zij to the point Zi'^j> immediately on its left, on its right, above 
it, below it, respectively. The first idea was to define: 

= TTi ^ + 

^ ' r + l u + d 

2 2 
Since the solutions of A(j) = are the same if we multiply A by any function, I 
preferred to multiply this expression by ^ x to get an expression that, you 
will notice, does not change when we compose the local coordinate with a C-afHne 
map. So let: 

u+d u+d 
[Aujij = —^(Ui+ij - Uij) ^y^hj ~ Ui-ij) 

I + r I + r 



in dK, which is called the Harmonic measure with respect to 00. It is also the unique non 
atomic measure with mass 1 and support in K (in fact dK) that minimizes the energy Eim) = 
— JjfyK log \z — w\dm{z)dm{w) . The minimal energy is equal to — logp, and is called the capacity. 
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Inside Sq we get (ui-ij + Ui+ij — 2uij)/K + (wij-i + utj+i — 2uij) x K and 
outside Sq we just get the standard Laplacian. On its boundary it is a little bit 
more complicated. Now for corners: we take the standard Laplacian. It is not quite 
correct, but the heuristics is that getting the formula wrong on a bounded number 
of points will not influence too much the result, as long as we take something 
reasonable: here {Au)ij vanishes if and only if the central value is a barycenter of 
the neighborhing values, (the weights of the barycenter are allowed to depend on 




Figure 5: The result of the computer experiment. Left: a grid, which within the 
square Sq has been divided into rectangles of ratio 2. Right: the image of the left 
grid by (j). The rectangles have become conformal squares. Note: The grid G on which 

the PDE was solved to produce this picture is finer than the one which is drawn. The spacing of 
G is 1/64 times the length of the side of the square supporting the Beltrami form. And G extends 
beyond the image: it extends in the 4 directions to a distance of 10 times the square's side. 

Let me now describe the numerical method I used to solve the discrete equation 
Afjj = with the computer I used the Jacobi Relaxation Method: recall that Au 
vanishes at the vertex if and only if barycenter of the value of the 

neighboring vertices (with weights that depend only on The idea is then to 

replace Uij by this barycenter. More precisely, define inductively a sequence u^"^ of 

functions u^fj as follows: begin with uf'j = Zij, which agrees with imposed 

boundary values. Now given u^"'\ if (i,j) is not a boundary vertex, let w^"^^'' be 
the value of the aforementioned barycenter computed from the values of m^"' at the 
vertices neighboring («, j). If is a boundary vertex, define '"i"^^'' = "i"'*- The 
reader will notice that since we change the value everywhere (but on the boundary), 
the Laplacian of still not 0. However, the values of u^"-' are supposed to 

converge to the solution of Au = 0. 

Remark. I made no effort at trying to prove the the convergence of this scheme 
and the correctness of its limit, nor at estimating the speed of convergence. 

As this was a little bit slow, I first solved the equation on a coarse grid, and then 
used the (near) solution u^^\ for some big N, to define a starting point w^"^ for 
solving the equation on a finer grid (with the grid unit divided by 2, the value of w^'^^ 
on the newly introduced points being defined by averaging on its neighbors), and 



Acknowledging I am not an expert in numerical solutions of PDEs, there is probably room 
for improvement there. 
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went on refining the grid further several times. This resulted in a great acceleration 
of the convergence. 




Figure 6: The image of the square, as approximated by the program, for K ~ 2, 5, 
10, 20. 

Let us list a few defects of the method: 

• The definition of the modified Laplacian makes use of the fact that the side 
of the square and the major axis of the ellipses are parallel to the coordinate 
axes. It is not clear how to adapt it to more general ellipse fields. 

• It is computer intensive. 

• It may have a too high algorithmic complexity (the time it takes to reach 
a given accuracy e). 

• The grid is uniform, whereas a non uniform mesh would be more efficient. 

• When K tends to infinity one needs at the same time a finer grid and a 
bigger N. The Laplacian approach does not seem well suited to big values 
of X. 

The reader will find on figure [5] what the same program yields when we increase 
K. It was hard to get something convincing in a reasonable amount of time for 
K ~ 50. A few hints for improvement could be: 

• One could have used the Conjugate Gradient Method (see [J]) instead of 
Jacobi Relaxation. 

• One may use a better approximation than (t>{z) — z on the boundary of the 
big square, allowing to take a smaller one. 

• Replacing, within the square, the nxn grid by an nxm grid, with n/m « K . 

• Using a Finite Elements Method. 
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All this discussion was about solving Acf) — 0. Amongst other approaches, there 
are: 

• Circle Packings (see [H] and [B] for instance). However, it probably needs 
a very dense mesh too when K tends to infinity. 

• Fourier series and the Hilbert-Beurling singular integral operator. See for 
instance [D], [GK| . One could also periodize the problem (using the big 
square as a fundamental domain) and use standard 2D Fourier series on 
the torus. 

• Discrete Riemann surfaces. See [Ml. 

But I did not push my experiments further. In any case, big values of K seemed 
to require a lot of work and computer time. The behavior of (/>(Sq) when K — > -\-oo 
remained a mystery for me. I was short of theoretical arguments, and there was no 
hope of an explicit formula there, unlike the case of the disk. . . Or so I thought. 




Figure 7: Five guesses of the possible limit shape, if there is one, together with a 
possible way to tend to it. 

2.3. A few guesses. Imagining that 0(Sq) has a limit when K — > oo, what 
could it be? It should be noted that the spiral, which was quite slowly turning for 
K — 2, does it faster and faster. But does its diameter tend to 0? Figure [7] gives 
a few possibilities. We will see later in the article that the reality is even more. . . 
interesting. 

3. An explicit formula 

It turns out that there is an (almost completely) explicit formula for (f>^^. It 
allowed the author to draw very precisely on the computer the image of the square, 
and to push the value of K up to 10^° and beyond. . . Thanks to these pictures, 
it was possible to guess the limit of 0(Sq), at least qualitatively. The formula can 
be proved to have a limit that allows to give a guess of the actual limit of (/)(Sq). 
We will see that the theoretical framework used to discover the formula gives an 
interesting interpretation of these limits. 
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Remark. What we will explain here generalizes to the straightening of Beltrami 
forms defined on C, piecewise constant, with polygonal pieces, and finitely many of 
them. This will be covered in a forthcoming article. 




Figure 8: The affine surface behind the formula. Case K — 2. The arrows are 
C-afRne maps (similitudes). 

Defining an affine surface: Consider the following Riemann surface: one chart 
is given by the exterior of Sq, union 9Sq, minus the four corners. The second chart 
is a rectangle given by the image of Sq minus its corners hy x + iy i—f x + i{y/K). 
Now glue the four sides of the first chart to the four sides of the second chart by 
four affine maps, as illustrated on figure Not only we get a Riemann surface, 
but the transition maps are C-affine. So we get an affine surfac^ S. 

Completing the Riemann surface: As a topological surface S is homcomorphic 
to the sphere minus five points. They correspond to oo and the four corners of 
the square. Let us prove that we can extend the Riemann surface S to these five 
pointsU by providing adequate charts near them. The easiest one to add is oo: we 
introduce the chart B{0,e) and glue it to Chart 1 with the map z i-^ 1/z. To add 
a corner, remark that a neighborhood of it is isomorphic, as an affine surface, to a 
neighborhood of of the following affine surface: C with a straight slit radiating 
away from 0, the two sides of the slits being glued together by a map of the form 
z 1-^ Kz. We have already seen in section [TT] how to conformally map this to C*: 
a branch of the map z i—s- will work, with 

a = 1/(1 — \og{K) /2i7T) for the upper right and lower left corners, 

a = 1/(1 + log(/-i:)/2i7r) for the other two. 

This gives a conformal chart in which we can add the origin, i.e. the corner. There- 
fore we have proved: 

■^Normally, charts should be open. So what is meant is that for all sufficiently small open 
sets containing our two charts, gluing them with the previous affine maps (more precisely by the 
unique C-affine extensions), we get Riemann surfaces that are all "the same", i.e. isomorphic by 
the obvious map. 

^An affine surface is a topological surface with an atlas whose transition maps are all locally 
C-affine maps, i.e. of the form z az + 6 (being locally C-affine on an open sot, the transition 
map is affine on each connected component of its domain of definition). An affine surface is in 
particular a Riemann surface by the same atlas. 

^Consider a topological surface S and a Riemann surface structure on an open subset U. 
Assume there is a point x Si S that does not belong to U but has a neighborhood included in 
U U {x}. It is not always possible to extend the Riemann surface structure to (7 U {x}. 
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Lemma 1. Each corner in S has a neighborhood that can be conformally mapped 
to a neighborhood of the origin in C, on which any branch of z ^ z^/" is an affine 
chart (for the value of a defined above). 

Uniformizing to the Riemann sphere: We have thus extended the Riemann 
surface S into a compact Riemann surface S, with five more points. Since S is 
homeomorphic to the sphere, it is conformaUy equivalent to the Riemann sphere S. 
Let $ : 5 — > § be the unique isomorphism normahzed by the foUowing condition: 

(1) = z + + o(l) as z — > oo 

where is the expressiorfl of $ in Chart 1. The subset 5 of 5 is mapped by $ to 
the complex plane minus four points. By the symmetries of the surface, the chosen 
normalization, and uniqueness of <I>, these points must be of the form zi = x + iy, 
Z2 ~ ~x + iy, zz — —X — iy, Z4 = X — iy, for some a; > and y > 0. 

Connection with the straightening: We can use the affine surface S, and the 
uniformization $ of its completion, to recover the straightening </> of the Beltrami 
form that is our object of interest. Let ^ : C \ {1 + i, —1 + i, —1 — i, 1 — i} — > 5 
be defined as follows: 4* maps a point z outside the square Sq to the point of S 
which has coordinate z in Chart 1. It maps a point z — x + iy within the square to 
the point of coordinates x + i{y/K) in Chart 2. Then $ o vj/ is a homeomorphism 
from C ^ C, that satisfies $ o '^{z) = z + + o(l) as z — > oo. It is a real 
diffeomorphism on the complement of the boundary of the square, that straightens 
the Beltrami form fi. By quasiconformal erasability of lines, it is the straightening 
of ^: 

(2) (jy^^o^. 

Equivalence of an affine surface structure compatible with a Riemann surface 
structure, and the distortion derivative: The affine surface structure on S, i.e. the 
notion of affine chart, can be transported by $ to an affine surface structure on 
C\ {zi, Z2, Z3, Z4}. Indeed, it is given by some atlas on S, and by composing charts 
with we get an atlas on C \ {21,22,^3,^4}. Now consider any affine chart 

{U, /) on C \ {zi, Z2, 23, Z4}, where U is an open subset of C \ {zi, Z2, 23, 24} and 
f : U ^ C The map / is analytic, because $ is. Since it is also injective, /' does 
not vanish. Since the transition maps between affine charts are affine, the complex 
number ri{z) = f"(z)/ f'{z), called the non-linearity or the distortion derivative of 
/, does not depend on the choice of the chart ([/, /) near z. Indeed for any two 
affine charts fi : Ui ^ C and f2 ■ U2 ^ C, and any z G Ui O U2, there exists 
a,b G C with a 7^ such that /2 = afi + b holds near z. Reciprocally, from ri{z) 
it is possible, at least locally, to recover the affine charts: / = / exp (J 77), the 
two integration constants accounting for a and b. In other words the function 77 
characterizes the affine structure. 

We therefore have a well-defined holomorphic function 77 : C \ {zi, Z2, 23, 24} —^ 
C, satisfying 77(2) = f"{z)/f'{z) for all affine charts / near z, and completely 
characterizing the affine structure. By equation ([!]), since f — (j>i is an affine 
chart, we get that 77(2) — > as 2 — > 00. 

Behavior of the distortion derivative under a change of variable: To study the 
singularity at Zi of 77, we will make a convenient change of variable: by Lemma [T] 
there exists a conformal map ip from a neighborhood of Zi to a neighborhood of 
in C such that for any branch 5 of z z^/", the map / = g o 7/) is an affine chart of 
C \ {zi, Z2, Z3, Z4} near zf, the number a equals 1/(1 — log{K)/2i7T) for zi and Z3, 
and 1/(1 + \og(K)/2iTr) for Z2 and Z4. The quantity g"{z)/g'{z) that characterizes 
the affine structure in the new conformal chart is equal to (3/ z where /3 — — 1. 



'i.e. 01 = <& o chartj ^ where charti goes from the corresponding open subset of 5 to C \ Sq 
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At Zl and Z3, (3 = — \og{K)/2i7T. At Z2 and Z4, /3 ~ log(iir)/2i7r. Since f = g o 
we have the composition formula: 

f" a" lb" 

(3) 77 = ^'x^°^ + 77- 

Now ip' does not vanish at the origin, and it follows from this and from g" / g' = (3/ z 
that the map r] = f" / f has a simple pole at zt with residueQ (3. 

The distortion derivative is completely determined: Indeed 77 is a holoniorphic 
function over C \ {zi, Z2, Z3, 24} that tends to at 00, and has simple poles poles at 
-Zl) Z2, Z3, z^ with residue — log(i^)/2i7r at z\ and Z3 and log(X)/2i7r at 22 and z/^. 
Such a function is necessarily rational and: 

4 T]{z) = + 

2z7r \z — z\ z — Zl z — zz z — Zi ^ 

So we know 77 explicitly (and thus the affine charts /, and from this the inverse 
of the straightening) provided we can determine the value of Zi . This can be done 
numerically, far much faster than for a Partial Differential Equation, as we will see 
in section H.3. 41 

Local expression of the affine coordinates: Explicitly solving the equation 
/"//' = yields: 

log K / \ 
log/' = -rr. — log(z - Z4) - log(z - Z3) + log(z - Z2) - log(z - Zl) + const 
2i7r V / 

\0gK ( {Z - Z2){Z - Zi) \ , 

— ■ log — + const 



2iTT \ (z — Zl)(z — Z3) ^ 

where we can guess there will be branch problems. At least locally the charts (C/, /) 
on C \ {zi, Z2, Z3, Z4} have the expression 

log K 

(5) i )^'^ .z+. 

J \{z - zi){z - z^) J 

for some a G C* and 6 G C. Determining the two affine charts corresponding to 
Chart 1 and Chart 2 thus requires to determine the value of a and b for each. 

Going global: The statement above is local. To extend it a little bit, we must 
be careful about monodromy questions. For one thing, the integrand 

\{z~ zi)(z - Z3) 

is a multivaluated expression: the set of possible values at a given z is of the 
form {A'"v I n G Z} for some v £ C* . Also, any affine chart / with a connected 
domain of definition extends to any bigger and simply connected open subset of 
C \ {zi, Z2, Z3, Z4} into a unique solution of /"//' = but the extension is not 
necessarily injective. 

Extending an affine chart along a path: Consider a path that is contained 
in C \ {zi, Z2, Z3, Z4}, ends included. Let zq be its starting point. Consider the 
initial values /'(zo) and /(zq). We can choose any determination mo of m(zo) 
and let b — /(zq) and a — /'(zo)/too- Then the integrand m(z) can be followed 
continuously along the path and its integral along the path is well defined, which 
gives a certain value for the expression a J m(z)dz + b. Any homotopic path will 
yield the same value. Let us call this the continuation of the affine chart / along the 



^In fact the equality of the residues of /"//' ^^nd 9" /§' still holds if the pole has higher order 
or is an essential singularity, as can be seen by multiplying equation ^ by dz and integrating 
both sides on a small loop around z;. 
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path. It depends only on f{zo), /'(^o) and the homotopy class of the path@ If the 
path is completely contained in the domain of definition Def (/) of an afline chart / 
(or the domain of definition of an analytic continuation / of an afhne chart), then 
we recover the value of / at the end of the path. This holds even if this domain is 
not simply connected: the mere existence of / implies that there is no monodromy 
within Def(/). In particular, there exists a branch of m defined on Def(/). If the 
path is not completely contained by just homotopic to a path completely contained 
in Def(/), we recover the same value. 

Expression of Chart 1 : Identifying S with C \ {zi, • • • ,24} via $, Chart 1 can 
be considered as defined on Ui — (/)(C\Sq) and equal there to (f>~^, since vj/ = id on 

the complement of Sq. Let us choose the branch of z 1-^ m{z) — ^ [^I^^j[^Z^!^| ^ 

on Ui that tends to 1 at 00. According to the above analysis, (/)~^ = b + a x { 
a primitive of m on Ui). Since (p~^'{z) — > 1 at 00, we must have a = 1. Note 
that to(z) — 1 is necessarily integrable on pathes within U that tend to cx), because 
4>~^{z) — z has a limit at 00 and m{z) = (In fact, m{z) = 1 + 0{z~^) at 00, 

as can be shown by a direct Taylor expansion, using the fact that Z2 + Zi = and 
zi + Z3 = 0.) Since ^^^(z) — z tends to 0, we get: 

Proposition 2. The straightening (f> of the Beltrami form defined in section \2.1\ 
has the following property outside of the square Sq; for all z ^ Ui = 0(C \ Sq), we 
have 

r-M^. +/'(( '- i"-"! )"^-')^"- 

J 00 \ \{u ~ zi){u - z:i) J J 

on any path from 00 to z that is contained in Ui ( or just homotopic within § \ 
{zi, Z2, Z3, Z4} to a path contained in Ui with the same endpoints), for the branch 

of the expression I j defined on the path and equal to 1 at 00. 

Remark. There is a difficulty in using this statement: it does not tell us how to 
determine whether or not a given z belongs to Ui. Moreover, remember that Ui 
spirals around the points Zi. Therefore, given some z € Ui, it is not obvious how 
to determine if a path going to 00 is homotopic in the complement of the z^'s to 
one staying in Ui. 

Expression of Chart 2: Using Chart 2 translated by i{l — K^^), we can extend 
Chart 1 to f7 = (j){{C \ dSq) U (-1 + i, 1 + i)). It corresponds to gluing Chart 2 to 
Chart 1 along the upper side of the square. According to the analysis above, the 
branch of m we have chosen on Ui extends to U, since the affine chart does. Let 
01 be the unique argument of zi that belongs to (0,7r/2). Notice that at 2 = 0, 
the quantity (^-^2)(^-^4) ^^j^^g ^^^^^^ ^-i4e, ^ rj.^^ p^^j^ 

ffoine from oo down to 

\yZ z\)yz z^ ) 

within the imaginary axis is contained in U. Following the logarithm of 

along the path, starting with the determination of log 1, leads to the value —iWi. 

From this, it follows that m(0) = exp(-i46'i \og{K)/2iT:) = exp(-26'i \og{K)/'K). 

Proposition 3. For all z E U2 = (f>{Sq), we have 

r-M^^(/Y !°-ii"-"i )^^A 

\Jo \{u- zi){u- Z3) J J 

where A(x + iy) — x + iKy, and the path of integration between and z is contained 
in U2 (or just homotopic within §\ {zi, Z2, z^, Z4} to a path contained in U2 with 



Thus we have a well defined analytic map on the universal cover of C \ {21, 2:2, 23, 24}. It is 
called the developing map of /. 
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r 



Figure 9: Extending Chart 1 to = 0((C \ aSq) U (-1 + 1 + i)). 

same endpoints ), for the branch of the expression ( j defined on 

the path and taking at z = the value exp(— 20i log(_R')/7r), where 9i is the unique 
argument of zi that belongs to (0,7r/2). 

Remark. We have the same difficulties as in the previous proposition. 

In the next section, we will explain how to practically use these formulae and 
how to determine the value of zi and the set 0(Sq). 

4. Practical aspects 

4.1. Integrating towards a singularity. Consider the continuation of an afhne 
chart along a path, as defined in the previous section. The following lemma tells 
that under certain conditions, the path can tend to one of the corners Zi. 



Figure 10: The half plane Re (z) — rim (z) < 0. Left: when r < 0. Right: when 
T > 0. A lift L{t) of the path 7(t), satisfying exp(i(t)) = 7(t) — Zi is drawn. The 
real part of the lift of the path must tend to — oo and its distance to the line must 
tend to +00. 

Lemma 4. Consider an initial point point zq G C \ {zi, Z2, Z3, Z4}, and a path 
7 : [0, 1] — > C from zq to one of the corners Zi, that is contained in C\ {zi, Z2, Z3, Z4} 
except at its arrival. Let L : —t C be a continuous lift such that exp(L(t)) = 

7(t) — Zi. Let T = \og{K)/2'K if i ~ 1 or 3, and r = —\og(K)/2TT if i — 2 or 4. 
Then a continuation of an affine chart along 7 will converge at t — 1 if and only if 

Re(L(t)) - Tlm(L(t)) — > -00 
t^i- 

(see Figure]!^. 

Consider two paths starting from zq, ending at Zi, satisfying the above condition 
on their lifts, and homotopic within C \ {zi, Z2, Z3, Z4}. The continuation of an 
affine chart along these paths will converge and tend to the same value at t = 1. 
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Proof. We may post-compose the affine chart with an affine map. 

First claim: Let /3 = it. We saw that there is a change of variable z = il}{w) 
defined near z^, sending Zi to and such that, locally, the affine charts g — f ° 4' 
satisfy g" [w) / g' [w) — (ijw (and thus are are branches oi w ^ aw^^^ + fe). This 
affine structure can then be lifted by the exponential: let w = exp(u) and h — goexp. 
Then h" /h' — (3 + 1. This equation has global solutions: let us choose h{z) = 
exp + l)^). Let W{t) be a continuous lift such that exp(W^(i)) = %l)~^{'^{t)), 
defined for t close enough to 1. Then /(7(t)) = h{W{t)) = exp((/3 + l)W{t)). 
Since ij) is an analytic diffeomorphism, W{t) — L{t) + const +o(l) as t — > 1. 
Since | exp((/3 + 1)VK(0)| = exp (Re ((/? + l)VF(t))) = exp (Rc ((1 + ir)VF(t))) = 
exp (KeW{t) — rim the first claim follows. 

Second Claim: by homotopy, it is enough to prove it for paths contained in a 
neighborhood of z,;. In that case, we can apply the change of variable Tp^^ as above, 
from which it follows that both continuations of / tend to 0. □ 

For instance, a curve tending to some Zi and having a tangent there will satisfy 
the hypothesis of the lemma. 

4.2. The parameter problem. We will denote by log^ the principal branch of the 
logarithm: it is defined on C \ (— oo, 0] and takes values into the strip M x (— tt, tt). 
Recall that we have defined a multivaluated function 

(Z - Z2){Z - Z4) 



•m{z) 



[z - zi){z - Z3) 



4.2.1. Theory. If instead of the unit square Sq on which the Beltrami form is 
constant, we decided to take a rectangle R defined by —h < Re (z) < h and 
—V < Im (z) < V, we could still define an associated affine surface and confor- 
mally uniformize it on the complement of four symmetric points in C, with an 
isomorphism normalized at 00 by z S Chart 1 i— > z + + o(l). This would give 
another value of zi, apart from which the affine charts would have exactly the same 
formula, with the same exponent log{K)/2iTT. This allows to define a map 

Z(h + iv) ~ z\ 

that depends on K. 

Independently of any choice of a rectangle, consider a given complex number 
zi in the upper right quadrant. Consider the path 7 from 00 to z\ expressed by 
zi + (1 + i)r with r G (0, +00] going from +cx) to 0. Let 

S(zi) = z + (to(u) — l)du 

where m(z) denotes the continuous branch along 7 that equals 1 at 00. This map 
also depends on K. 

If zi = Z{h + iv), then by Lemma|4]we must have S(zi) = h + iv: indeed let 72 
be the image by (j) of the path h + iv + (l + i)r, with r g (0, +00] going from +00 
to 0; first the path 7 is indeed homotopic in C \ {zi, Z2, Z3, Z4} to 72, because both 
are contained in the upper right quadrant minus {zi}; second the condition on the 
lifts L{t) are satisfied: the lift of 7 is horizontal and the lift of 72 is asymptotic to 
a perpendicular to the line Re (z) — rim (z) = 0, as can be seen using Lemma [T] In 
other words 'EoZ = id on the quadrant. In particular, Z is injective, S is surjective. 
The map Z and S are easily seen to be homogeneous: for all z in the quadrant, 

VA > 0, Z(Az) = XZ{z) and S(Az) = AS(z). 

Note that Z{h + iv) is the image of the upper right corner of the rectangle R by the 
straightening. It is thus a continuous function oi h + iv, and it extends continuously 
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on the boundary of the quadrant to the identity (because it then corresponds to a 
null Beltrami form). By this and homogeneity, Z must be surjective, and tend to 
cxD at oo. The closed quadrant union oo is a compact set on which Z is bijective, 
continuous and has an inverse S. This implies that S is continuous. Therefore, Z 
and S are inverse homeomorphisms. Let us sum-up: 

Proposition 5. The function 'E. is a homeomorphism of the right quadrant, with a 
continuous extension to its boundary equal to the identity there. 

4.2.2. A method. Here we present a numerical method to obtain zi. There are 
probably many other ways of doing this. 

To determine zi, and thus the other Zi by symmetry, we want to solve for 

S(zi) ^l + i. 

Let us choose a branch of the integrand function m(z): let 



mp(z) — exp 



{~^~~i}^^p^^' - Z2) + logp(z - Zi) - logp(z - Zl) - logp(z - Z3)) 



where log^ denotes the principal branch of the logarithm. It is naturally defined 
on the complement in C of the two horizontal half-lines on the left starting from zi 
and Z4: 

Dcf (mp) = C \ ((zi + (-C50, 0]) U (Z4 + (-00, 0])) . 
See figure [TlJ Near 00, it tends to 1 and thus coincides with the branch used 





k 

Zl 








► 


2.3 


Z4 



Figure 11: The domain of definition of mp 



in Proposition [21 which is equal to 4>~^' . Let us consider the continuation fp to 
Def(mp) of the afhne chart near 00: fp{z) = z + J^{mp{u) — l)du for any path 
of integration contained in Def(mp) U {00}. Thus 

(6) /p(z') - fp{z) = mp{u)du 

for any path of integration from z to z' contained in Def(mp). Note also that any 
path tending to zi within Def(mp) will satisfy the hypotheses of Lemma|3]and thus 
the formula defining fp{z) also works for z = zi, and we have moreover 

/p(zi) = S(zi). 

To better avoid numerical instabilities, especially when logK gets big, instead of 
computing /p(zi) directly, in the author ended-up after a few attempts with the 
following way of determining /p(zi). There is no claim that this choice is optimal or 
even good, but it worked in the numerical applications. Let zi — x+iy, Z2 — ^x+iy, 
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Z3 = —X — iy and Z4 = x — iy. Compute 



lb 

Ic 



Ic 



a 



.fp{zi +i) - fpiiy + i), 
fp{zi) - fp{zi 
fp{zi+i)- fp{x + l), 



each of these numbers being determined by numerical integration of equation ^ 
on the segment between the two points where fp is evaluated. See figure \Vl[ Then 



So we have a method to compute S. Recall that this function is homogeneous. A 
method to find the correct value of zi, the one for which = 1 + i, is then the 

following: 

• find the unique 9i £ (0, tt/2) such that the complex number c — S(exp(i6'i)) 
is collinear to 1 + for this we may look at the slope s = Re (c)/Im (c) and 
use the secant method to solve equation s = 1; 

• set Zl = exp(j6'i)/Rec. 

4.2.3. Examples of parameter determinations. To compute the integrals, the Simp- 
son rule has been used, with an adaptive subinterval size, except near the singularity 
Zl, on which the integrand has been approximated by a constant times a complex 
power of the remaining distance, which can be integrated explicitly. Let us be more 
precise. 

Recall that we have chosen a branch of the integrand function m(z): 



mp{z) = exp — — ( log (2; - Z2) + log (z - 24) - log„(z - zi) - log (z - Z3) 



where logp denotes the principal branch of the logarithm. 

The adaptive method consists in subdividing an integration segment into unequal 
subintervals. The subintervals are constructed one after the other, in increasing 
order on the segment. The size of each subinterval is decided once its starting point 
is known. The error in Simpson's method, which consists in estimating / by 
(/(a) + 4/((a + b)/2) + f{b)) /6, is bounded from above by (b - a)^ sup |/2880. 

A very crude estimate was used for sup Iwp''^ (z)| on [a,b]: M\mp{a)/\a — zi|^ for 
some constant M. The integral on [a, b] being close to (6 — a)\mp(a)\, the relative 



Figure 12: The integration segments. 



Re {la + h) + ilm [h + h) = 
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Table 1. 



K 


Rc(zi) 


Im (zi) 


do 


To 


steps in J 


2 


1.2480 


0.7676 


0.2 


0.01 


38 


10 


1.6483 


0.4276 


0.2 


0.01 


37 


100 


1.8378 


0.2350 


0.2 


0.005 


39 


10^ 


1.8976 


0.1184 


0.2 


0.002 


39 


10^ 


1.9066 


0.0789 


0.2 


0.001 


46 


10^ 


1.9104 


0.0526 


0.2 


5.10-4 


49 


1012 


1.9117 


0.0394 


0.2 


5.10-4 


49 


102" 


1.9127 


0.0236 


0.2 


5.10-4 


49 


10™ 


1.9132 


0.00947 


0.2 


2.10-4 


53 



error is then estimated to M 

a—zi 

is small. This amounts to taking 



4 



/2880. Then b is chosen so that this quantity 



b — a — do X |a — zi| 

for some do > 0. 

For the segment that goes to zi, we fix some distance tq > to zi below which 
we stop using Simpson's rule, but instead do the following. Assume that we want 
to integrate mp{z)dz on the segment from zi + w to v. Depending on v, there is a 
constant mo € C* such that mp(zi + tv) ~ mo exp(— /31ogt)) as t — > 0+, where 
P = \ogK/2iTr. The right hand side of this equivalence is easy to integrate on the 
segment from zi + v to z\: J mo exp(— /31og(i)))dz = vmo exp{—(3\og{t))dt = 

—vmo/{l — 13). In our case, v — roi- Thus nio = exp (^f3(^\ogp{zi — Z2) + logp(zi — 

Zi) - logp(z3 - Z2) - logp(roi)) j . 

In Table [T] we show the results of some computer experiments. The values of 
do and ro have been chosen by hand: the error has not been estimated, even if 
it is likely that the values are correct up to adding ±1 to the last decimal. By 
convention all decimal numbers are approximations by truncation (=by default). 
The field "steps in J" refers to the total number of subintervals the three integration 
intervals have been subdivided into. It does not count the subinterval that touches 
Zi . The values of zi and "steps" shown in the table are those corresponding to the 
last iteration of the secant method. The secant method never required more than 
20 iterations, usually less (depending on the initial guess, obviously). 

As expected, Im(zi) — > Oa.sK increases. It does so quite slowly: it seems to 
be comparable to l/\ogK. It is not completely surprising as we know that in the 
formula for m the constant K appears only as logi^T. As for Re(zi) it seems to 
tend to a value close to 2, but < 2. 

In these computations, the sufficient number of integration steps to reach a 
precision of 10-4 is surprisingly low. The same experiments have been run with 
lower values of do and ro, whence more integration steps, and yielded the same 
10~4-approximations: on Table [21 we see that the value stabilizes, and that 11 
decimal places could already be obtained with only approximately 1500 integration 
subdivisions. Similarly, for the extremely big value K — 10^", Table [3] shows that 
less than 2100 steps were enough to give an approximation with an error < 10~^^. 
Again, in the absence of a abound on the error term, these experiments are given 
for information, not proof. 
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Table 2. Numerical experiments for K = 2. 



Re(zi) 


Tm ( z^ ) 


dr\ 


' U 


steps in J 


1 25284529816 


764992357071 


1 


0.1 


5 


1.33602299688 


0.722976075222 


1 


0.001 


5 


1.15536861038 


0.460939047597 


1 


10-12 


5 


1 23452767620 


766441199359 


0.5 


0.5 


9 


1 24582003447 


767669175930 


0.5 


0.2 


10 


1 24756645844 


767688650200 


0.5 


0.1 


11 


1.24815951395 


0.767643162006 


0.5 


0.01 


14 


1.24816586438 


0.767640989544 


0.5 


10-4 


21 


1.24816586643 


0.767640982471 


0.5 


10-8 


34 


1.23447420199 


766433209762 


0.2 


0.5 


21 


1.24748486405 


0.767688303419 


0.2 


0.1 


28 


1.24807062425 


0.767646071328 


0.2 


0.01 


38 


1.24807622667 


0.767644491781 


0.2 


10-4 


59 


1.24807()227()5 


().767()44491379 


0.2 




100 


1.24748382611 


0.767688184430 


0.01 


0.1 


564 


1.24806951508 


0.767645985170 


0.01 


0.01 


793 


1.24807511121 


0.767644410897 


0.01 


10-4 


1251 


1.24807511157 


0.767644410561 


0.01 


10-8 


2167 


1.24748382610 


0.767688184430 


0.001 


0.1 


5633 


1.24806951508 


0.767645985169 


0.001 


0.01 


7934 


1.24807511120 


0.767644410897 


0.001 


10-4 


12537 


1.24807511157 


0.767644410560 


0.001 


10-8 


21743 


1.24807511157 


0.767644410560 


10"^ 


10-8 


2 175 079 


1.24807511157 


0.767644410565 


0.01 


lo-'' 


1480 


Table 3. Numerical experiments for K = 


iriSO 


Re(2i) 


Im(2;i) 


do 


To 


<^tpTT^ in / 


1.91325679285 


0.00947278823076 


0.2 


2.10-4 


53 


1.91325221530 


0.00947281249592 


0.2 


10-8 


Q7 


1.91325221554 


0.00947281249722 


0.2 


10-12 


138 


1.91325400121 


0.00947285755832 


0.1 


10-8 


202 


1.91325400121 


0.00947285755831 


0.1 


10-12 


290 


1.91325406829 


0.00947285956006 


0.05 


10-8 


414 


1.91325407082 


0.00947285966708 


0.02 


10-8 


1045 


1.91325407086 


0.00947285966944 


0.01 


10-8 


2098 


1.91325407086 


0.00947285966959 


10-4 


10-12 


302810 



4.3. Using the formula. Recall that we have defined a multivaluated function 

logK 

{Z - Z2){Z - Z4) \ 

{z - zi){z - Zs) J 

4.3.1. Computing the inverse of the straightening. Now, let us describe one way to 

determine whether a given point z belongs to (/)(Sq), and what is the value of 4>^^(z). 
Consider any path 7 from 00 to z. Follow along this path the branch of m(z) that 
equals 1 at 00. In order to follow (j)~^{z), begin by integrating m{z) — 1 along the 
path. Then (j)~^{z) = z + J^(m— 1) at the beginning, i.e. as long as (j)~^{'-f{t)) ^ Sq. 
The fact that the path starts from 00 is not a problem for the integration, because 
the integral is convergent near 00 (if one prefers, applying the change of variable 
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z ^ 1/ z leads to a standard path integral). Another possibility is to determine the 
value of (t)~^{zQ) for one point zq € C\i?(0, 2), in which case zq ^ 0(Sq), using either 
a path starting from oo, or a power series expansion of the formula, and then use zq 
as a starting point for the path 7: 0~^(7(t)) = 0^^(zo) + /J^*^ m{z)dz as long as the 
computed value of 4'~^{j{t)) does not belong to Sq. Whenever the computed value 
of (j)~^{'y{t)) enters in Sq, the integration goes on but with an integrand equal to 
A{m{z)dz) instead oi'm{z)dz, where A{x + iy) = x + iKy. Whenever the computed 
value of (j)~^{7{t)) gets out of Sq, we switch back to the normal integrand. And so 
on, up to the end of ^{t). We could have started at instead of 00, since we know 
that ^(0) = 0. In this case, begin by integrating A{m{z)dz), and follow the branch 
of m{z) that has the correct initial value (which is given in Proposition [3]) . 

4.3.2. Computing the the straightening. Computing (^{z) looks a lot like solving 
an Ordinary Differential Equation. Recall that cj) = (f> o \I>. Compared with the 
computation of 4>~^, the difhcult problem of determining whether a point z G C 
belongs to </'(Sq) is replaced by the trivial problem of determining whether z € Sq. 
So there is two cases: either z ^ Sq, in which case ^'(z) is the point with coordinates 
w = z in Chart 1. Or z G Sq in which case ^'(z) is the point with coordinates 
w = Re (z) + ilm {z)/K in Chart 2. 

Choose a point wo in the same chart for which zq = <&(wo) and <i>'(wo) are already 
known, choose any determination of m(zo) and choose a constant a e C* such that 
l/<i>'(wo) — am(zo). For instanceQ in Chart 1 we may choose wq = 00, for which 
Zq = $(wo) = 00 and ^'{wq) = 1, and choose m(zo) = 1 whence a = 1. In Chart 2 
we may choose wq = 0, for which $(wo) — and $'(wo) = exp(20i log{K)/n) with 
the notations of Proposition [3 and choose m(zo) = exp(— «46'i log(ii')/2i7r) whence 
a = 1. 

Choose a piecewise linear path 7 from wq to w that does not leave the chart. 
Then 



dt am($(7(i))) 

where the correct determination of m at $(7(t)) is obtained, as usual, by continu- 
ously following a branch along the part of the path ^{'-/(t)) under construction. 

4.3.3. Drawing the boundary of the square. The previous method can be used to 
draw the boundary of the square: first determine 0(1) and 4>{i), then draw the 
image of the path from i to 1 + i, do the same for the path from 1 to 1 + i, and use 
the symmetries to draw all the boundary of the square from these two curves. 

However, this method is a bit unstable for drawing (f){[i, ! + «]) when K takes very 
big values. The reason is that the distance between (j){i) and (t>{0) is comparable 
to 1/K. (Also since 0(1) is only approximated, it means that the solution of the 
corresponding differential equation will miss zi, even if slightly; since for big values 
of K there is a lot of spiraling around zi, could this yield big inaccuracies?) The 
author ended up following the two paths in the opposite direction, i.e. from 1 + i 
to i and from 1 + i to 1. We cannot begin right on 1 + i, because m degenerates at 
zi = 4>{i + 1). So we need to begin slightly off: the difhculty is then to determine 
the initial values of $(z) and <I>'(z). 

Remark. It is not really important to start right on, or extremely close to, the 
spiraling image of the side of the square by (p. Indeed, we are not at a fixed point of 
the vector field, but at another sort of singularity. Here small error is expected to 
remain small enough, see figure fT51 The problem comes from choosing the correct 



^Another possibility is to let wq = "J" ^(^o), for some chosen zg, where <t> ^ and m{zo) may 
be computed with the method of section |4]3TT] 
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Figure 13: The field lines of dz/dt = for a purely imaginary /3 are well defined 
as non-parameterized curves. They are shown for a small value of (3 on the left, big 
on the right. The field lines passing close to the singularity stay close to the one 
passing through. 



branch of m{z). The fact that in our case /3 is purely imaginary will make life 
easier. 

For this, let us use the fact that (f>~^ is, along the curve (p^i, 1 + i)), a primitive 
of a branch of m{z). To make this branch more explicit, let U = the upper right 
quadrant minus {zi}. The restriction of m to f/ is multivaluated, and is well defined 
on the universal cover U of U. There, it satisfies 

m{z) = ^{z)ex.p {(3\og{z - zi)) 

where /? = —\og(K)/2iTr and where ^ is a holomorphic function defined on the 
upper right quadrant, and 

= exp (^-^:^(logp(^i - ^2) + logp(2i - Z4) - logp{zi - 23)) 

where log^ denotes the principal branch of the logarithm. A primitive P of m is a 
multivaluated function on U, well defined on U, that satisfies 

P{z) = C{z) exp ((1 + /3) log(2: - zi)) + const 

where C is a holomorphic function of z defined on the upper right quadrant and 

co:=C(^i)=^(^i)/(l+/3). 

Take const = 0: then a point z E U on the curve (t>{[i,i + 1)) satisfies P{z) < 0. It 
turns out that, because /3 is purely imaginary, the condition P{z) < is independent 
of the lift of z G U inU, since changing the lift only multiplies exp ((1 + (3) \og{z — 
zi)) by > for some n G Z. From this it follows that the curve 4>{[i,i + 1)) 
spirals down to Zi asymptotically on the logarithmic spiral defined hy z = Zi+ re*^ 
and 

logii" 

TT = arg(co) +0+ logp(r) (mod 27r). 

One then starts from a point on this spiral for r small enough and follows the 
differential equation 

dz _ -1 

dt m{z) 

for some branch of m, until the imaginary axis is crossed. 



20 



ARNAUD CHERITAT 



Remark. The time ti — to taken to reach this axis will not necessarily be 1 but 
rather with n G Z, depending on which branch of m is chosen. The author did 
not manage to find a simple way to determine, given the initial point zq = z\ +re'^, 
which branch of to(zo) to take in order to get ii — to = 1- However, it is likely that 
there is one. 

The same can be done for the curve 1 + «)), with 

-| = arg(co) + 0+i^logp(r) (mod 2^). 

and 

dz —i 
dt m(z) 

4.3.4. Numerical experiments: a square on a diet. We present on figures [TJ] and [TS] 
the image of the boundary of the square under the straightening (jj, computed 
using the method of section [4.3.31 The differential equations dz/dt = —l/m{z) 
and dz/dt = —i/m{z) have been numerically solved with the fourth order Runge- 
Kutta method RK4, and an adaptive elementary interval size dt: more precisely 
\dt\ = eo\m(z)\, so that \dz\ is of order Eq. The value Eq = 10~^ was taken for the 
pictures. The list of points obtained approximates very well the curve, but is very 
big. A sublist was taken, and the corresponding polyline was drawn. 
What we observe, as K — > +oo, is that: 

• Re (zi) — > .Too > 

• Im(zi) — y 

• There seems to be a limit shape. 

• The limit shape is not a segment. 

• It has non-empty interior. 

• Yet it is not the closure of its interior: there is also a horizontal segment. 

• The limit shape separates the plane into several connected components. 

• The limit in the sense of Caratheodory of (/)(Sq) (which is the union of 
the limit shape and all the bounded components of its complement), is 
equal to [ 

— '^005 ^cxd] U U ( — Z?) where is a topological disk (a Jordan 
curve and its bounded component) contained in "Re (z) > 0" and with 
D n [0,Xoo] = {xoo}- It is smooth, except at Xoc where D has a cusp. Let 
us call D the bulb. 

• Near the endpoints of the segment, the boundary of the limit shape looks 
like a sheaf of tangent circles. 

• Near these endpoints, when the parameter K is big, the picture is reminis- 
cent of parabolic bifurcation. 

5. The limit 

5.1. Behavior of the distortion derivative. Let us note zi = x + iy with a; > 
and y > 0, so that 

Z2 = 0(— 1 + i) — —X + iy, zi = 0(1 + i) = X + iy, 
Z3 = (/'(—I ~ i) — —X — iy, Z4 = 0(1 — i) = X — iy. 

These numbers x and y depend on K. The previous sections suggest that, as 
K — > +00, y tends to 0, x tends to a positive value and has a non-affine limit 
on the complement of the square. 

Whatever value of K was chosen, the function is univalent on the complement 
of the square, thus in particular on C \ i?(0, v^), and normalized so that <f)[z) — 
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K=l 




Q D 




^ © 

Figure 14: The image of the square under the Beltrami form straightening, for 
different values of K. 

z + + o{l) as z — > oo, thus in particular (/)' — > 1 at oo. By a well known corollary 
of the area theorem 0(C \ B{0, \/2)) contains C \ B{0, 2\/2). Therefore first 

\zi\ < 2^2; 

second / = is univalent on C \ -8(0, 2a/2) and normalized like (f>: (f>^^{z) = z + 
+ o(l) as z — > oo. For such normalized univalent functions / there exists uniform 
bounds on the derivatives of a given order: in particular there exist functions Ao{R), 
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1 K = 10^°, closer view 



Figure 15: Zoom on the bulb. 



Ai{R) and A2{R) independent of / such that 

\z\ = R =^ Ao{R) < \f'{z)\ < A,{R) and |/"(z)| < A^iR). 
In particular \f"{z)/f{z)\ < A2{R)/Ao{R). Now, recall that 

1111 



+ 



Z3 



f'{z) 2i7V \Z — Z2 Z — Zi Z — Z\ 

Thus for all fixed 2; e C with \z\> 2^/2 one must have 

1 1 1 1 \ 
+ 0. 

^Z — Z2 Z — 24 Z — Z\ Z — Z^ J 

Since the Zi are bounded, this implies that collisions must occur: for all sequence 



+00 such that the four points Zi converge, we have 



2— lim Z2 



+ 



z— lim 2:4 
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z-iimzi ~ z-iimz3 ~ ^ whenever \z\ > 2^/2. Since the functions l/{z — a) for 
different vahies of a S C are hnearly independent on any infinite subset of C, then 
the foUowing unordered pairs are equah {hmzijhmza} = {hm Z2, Km Z4}. Using 
the symmetries between the Zi, we conclude that 

either x = Re (zi) — > or y = Im (zi) — > or both 

as K — > +00. We can even be more precise: since 

1 1 1 1 \ —Aixyz 



z-Zi z-zi z-z-i) ULiiz-Zi 



we have 



xy = o(- 



++00 WogK / 

However, it seems hard to prove with this method what is suggested by the computer 
experiments of the previous sections: that y tends to but not x, which shall tend 
to a positive value a; 00, and that has a non-afhne limit, thus rj has a non-zero 
limit, i.e. xylogK has a non zero limit. This is the object of the next section. 

5.2. Behavior of the image of the corners when the dilatation ratio tends 

to +00. The objective of this section is summed- up in Theorem |9l 

5.2.1. Motivating the estimation of some integrals. Note that if the claims made 
above hold, then the function 77 tends to 

lim(ylogi^) / -1 1 



This limit function defines an afhne affine surface structure over C\ {—Xoo, Xoo\, by 
expressing its distortion derivative. We may expect it to be a limit in some sense 
of the affine surfaces for K < +00. We may also look for a description of this new 
affine surface, in terms of charts and gluings. From the pictures, it seems that the 
part corresponding to C \ U {—D), where D is the bulb, is obtained by taking the 
complement of the square and gluing the upper and lower sides by a translation. 

Noteworthily, the function m corresponding to the limit function 77 would have 
no monodromy anymore: 

m(z) = exp 



(z ^oc) (-^ ^" '^oo) 

with s = lim(y log K) /tt. A primitive of m on the complement of [— Xqo, 2^00] U Z? U 
(— D) would be a conformal map to the complement of the square. Again from the 
pictures, and by analogy to the case K < +00 we expect that the pair (xoojs) is 
the unique pair of positive reals such that 



and 



height = [ {m{it) - l)d{it) = 1 

Jt=+oc 

width — I m{t)dt — 1. 



The rigorous justification will be a bit complicated. We will begin with a slightly 
different version of the second integral, for K < +00, and let K — > +00. 
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5.2.2. Estimating integrals. To estimate the limit of zi as K — > +00, let us con- 
sider the continuation along the path from +00 to within (0, +00) of the afSne 
chart near infinity. Because of the symmetries, the image by of the path 
ends up at the center of the square in Chart 2 (recall that $ conformally maps the 
abstract afSne surface to C \ {zi, Z2, Z3, Z4}). The map 

Charti o $-1 on C \ </.(Sq) 

{1- K) + K Charta o $-1 on (j){Sq U{i-l,i + 1)) 

satisfies the equation /"//' = r] and coincides with (/>~^ out of the square. It is 
thus equal on (0, +00) to the afiine chart continuation mentioned above. Therefore 

1-K = 0+ I {m{t) - l)dt 
J+00 

for the branch of m on (0, +00] that equals 1 at 00. Let us make this more explicit: 
with 

l{t) = logp(t - Z2) + logp(i - Z4) - logp(t - zi) - logp{t - Z3) 

where log^ denotes the principal branch of the logarithm. Using the symmetries 
between the Zi, this simplifies to 



l{t) = 1(2 argp(t -Z2)-2 argp(t - ^1)) 



where arg^ e (— tt, tt). And since t > and argp(t — Z2) € (—■ 7r/2,0) and argp(t — 
zi) € (-77,0): 

l{t) = i2e{t) with 0{t) = argp = argj, 

t — Zi Zi — t 

Recall that 

zi = x + iy. 

Then 

0{t) = argp -^-* + 'y = _ argp(f2 + y2 _ ^2 ^ -^^y-^ g (q, tt). 
X Xi ~\~ ty 

Thus 

e{t) = cotan-iX(t) with 

xit) = A + ify_^) 

2xy 2\x y' 

where cotan"^ takes its values in (0,7r). 

So 

K-l = (^e^p (^^^9{t)^ ~l^dt = j\---)dt + J^°°{---)dt. 

h h 

Bounding I2: We will use the following inequality: 



X>0 ^ cotan-i(X) < 



■X/] 

For l2,t> X thus X{t) > y/2x > 0, thus 
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Let US use the following convexity inequality: 

- 1 

0<s<S =^ -Ks- 



s ■ 

Applying this to s = ^-^/{l + X(i)/f ) and 5 = log(i4:)/2 gives 
Thus 



VK-1 

h < / —2 dt = 

1+- — +^-- 
TT \xy X y / 



r ^ dt 



l + - + ---) 1 + - — + — + - 

TT \ xy X y / TT \xy y X 



I VK — 1 , , , — s irJ-Kxy 



TT \xy X 

because J^°^ dx/{ax^ + b) = ^/Vab. Now recall that \zi\ < 2-^2, whence x < 2^/2 
and y < 2 a/2. Therefore 

h < {VK- 1)ttV8^. 

Bounding Ii : It is a lot easier: The integrand in /i is < iiT— 1 and the integration 
interval is (0, x) thus 

h < {K-l)x. 
Now recall that h + I2 = K - 1. Therefore 

liminf a; > 1. 

Let us use a slightly finer bound on the integrand: for t G (0, +00) we have X{t) S 
(X(0), +00) with X(0) = (I - |)/2. Thus 

h < {RP - 1)2; < (RP - 1)2^2 

withp = i cotan"^ -''^(0) < 1- This implies that p — > 1, i.e. that cotan~^ -'^(O) — > 
TT, i.e. that X{0) — > —00, i.e. that y/x — > 0. Therefore 

y — > 0. 

Making a pause: So far, we have proved that < liminf x and limy = 0. 
Recall that \x\ < 2\/2 and that a^ylogii' is bounded from above. We still need to 
prove that x has a Hmit and that xylogR has a positive limit. Let us pass to a 
subsequence Rn — *■ +00 such that x has a limit x ^-^id xylogR has a limit k. 
Necessarily x > 0- 

Proposition 6. 

1 



with a — , and 







X / exp ( -^—^ \du^\ 



X / ( 1 - cxp ( \\du=\ 



Proof of Proposition 
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Lemma 7. 



Proof. Wc omit the index n. 

K K 

with 



exp ^ — ^ — (^'(■") — tt)^ c?M 



2y 2'-x y' 

Now the integrand belongs to (0, 1) thus we may apply dominated convergence, 
and focus on the pointwise limit of the integrand for n — > +cxi. Since y — *■ and 
X is bounded from below, y/x — *■ and x/y — > +oo thus for a fixed u € (0, 1), 
§ + Mi - f ) - — > whence e{u) |^ + o{y/x) whence 

Since Kn — 1 = Ii + 12 and I2 = 0{\/Kn), Lemma [7] implies the first equality of 
Proposition [HI 

Let Y denote the positive part of the imaginary axis: Y = {iy | y > O}. Let 
I3 — J(m(z) — l)dz integrated along the path from 00 to within Y, using the 
branch of m defined on Y and tending to 1 at 00. Then by similar arguments as 
for Ii + I2, we have 

and 



idt 



/3 = ^ (1 - exp arg^ {x^ - y' + f ~ ^2xy) 

where argp denotes the principal argument € (— 7r,7r). 
Lemma 8. 

Proof. Let K = K,,. 

Is^i J ^1 - exp - ]ogK_ ^^^^^-1 ^ 

that we will write I3 — i /q^"" fn{t)dt, with 

zxy z V y a; / 

As t varies from to +00, the function X{t) varies from ^(^ — f ) to +00. For n big 
enough, — f ) > 0. Then, using 1 — exp(— w) < v for 1; > 0, and cotanw < 1/v 
for ^; > 0, we get that for n big enough, Wt > 0, < f{t) < f^^'^ll}"!^)^ ■ Since as 
n — > +00, xylogK converges and — — > > Oi the sequence of integrals 
I3 — i /g^°° fn{t)dt is dominated. So we may focus on the pointwise limit of fn{t)- 
The rest of the proof is similar to that of Lemma [71 with the change of variable 
t = xu and using cotan"-'^ f ~ l/ii as u — > +00. □ 

This gives the proof of the second inequality in Proposition [Sj 

End of the proof of Proposition 
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Let US recall the content of this proposition: for any sequence Kn — > +00 such 
that X has a limit x ^^id xy log K has a limit k we have 



X I exp ( -^—^ ]du^l 



with a — -2%, and 



Now for a fixed x the first integral is a decreasing function of a > and the 
second one is increasing. By dominated convergence the first integral tends to 1 
as fj — > and to as tr — > +00; the second tends to as cr — > and to 
+00 as a — > +00; both are continuous functions of a. Therefore their quotient 
second/first, independent of x, is an increasing continuous bijection from (0, +00) 
to itself. There is thus a unique value of cr > so that the quotient equals 1. There 
is thus a unique pair (x, tr) so that the pair of equations hold. There is thus a 
unique limit to x and to xylogK: 

Theorem 9. 1. There exists a unique solution {xoo,cf) G (0, +cxd)^ of the system 

W{Xoo,<7) = 1 



H{xoo,o-) = 1 



with 



W{xoa,<j) = Xoo I exp ( Y^-^ ) rfu, 





H{xoo,cr) = Xoo y ^1 - exp 
2. As K — > +00, 

TTtT 

zi > Xoo, Im(zi) 



77(2) 



2 log if' 
-1 1 



(z-.Too)^ (2 + Xoo)' 



Remark. If instead of the square Sq we had started from a rectangle of horizontal 
size 2w and vertical size 2h, then the same conclusion holds with the system replaced 
by 

VK(Xoo,ct) = W 

H{xoo,o-) = h 

5.2.3. Numerical values. Finding a and Xoo is similar to but simpler than finding 
zi. We obtain: 

limRe(zi) = Xoo = 1.9133480795... 

cr = 0.3628700796 . . . 

(JX 

— — = 0.3471483850... 
2 

lim (Im (zi) log A') = — ^ = 1.0905988161 . . . 

For instance, if we take K — lO^*', and the value computed for zi in tableOthen 
we obtain Re (zi) = 1.9132540708 . . . and Im (zi) log if = 1.0906032731 . . . which 
are close to the values above. 
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5.2.4. A remark about a residue. The following integral appearing in the second 
equation of theorem [5] can be expressed as a residue: 

H = H{xoo,cr) = Xoc J ^1 - exp _|_^2 j j 

Recall that at the limit the function m is not anymore multivaluated: 

s s 



m{z) — exp 



^ (-2' 3^oc) ^" ^oo) ^ 

with s = lim(y logisr)/7r ~ axo/2. Recall that ii? — J {m{z) — l)dz, integrated along 
the path from +oo to in the imaginary axis. By symmetry 2iH = J{m{z) — l)dz 
along the full imaginary axis followed downwards. Consider the integral J{m(z) — 
l)dz along the following curve: it starts from iR and descends the imaginary axis 
down to —iR; then it follows the half-circle from —iR back to iR and located on 
the right of the imaginary axis. This integral will be equal to 2i7r times the residue 
of m at on one hand. On the other hand, the vertical part tends to 2iH as 
R — > +00, and the circular part is 0{1/R) because m{z) = 1 + 0{z'^) as z — > oo. 
Therefore 

H ~ IT X residue(m at a;oo)- 
However, this residue is not really explicit. It can be computed using power series 
expansion and gives 

+ C30 



residue — > 

9 ^ \ A 



{2n)\ 



2 ^ V 4 / (7i!)2(n+ 1)! 

n=0 ^ y \ I \ ' I 



which converges quite rapidly (like the power series of exp(z)). 

But in fact, the same expression could have been computed directly from the 

definition of H as an integral: by expanding 1 - exp {^j^^ = - J2n>i ^ (tT^) ' 
integrating from to oo, permuting the sum and the integral (it is possible because 
the total mass is bounded), and shifting the index, we get 

+25 (-a)» du 

H = XooO ^ 



and the constant /q^°° ^-^^^^^„ can be computed explicitly, by induction (or as a 
residue. . . ). 

The author ignores if an analog approach can be set up for the integral W{xoo, c). 
5.3. The limit afflne surface. 

5.3.1. Introduction. Recall that any analytic function h on an open subset U of 
C defines an afhne surface structure on U, whose charts are the local solutions of 
/"//' = h (in particular they are analytic). We called h the distortion derivative 
of the affine structure. 

The function 77, which depends on K, is defined as the distortion derivative of 
some afRne structure on C \ {zi, 2:2, Z3, Z4}, which was obtained by conformally 
uniformizing an affine surface S defined by polygonal pieces in C glued along their 
edges by affine maps. It is thus natural to interpret its limit as K — > +00: 

axr^ / -1 1 

Voc(Z) 



2 \iz - Zoc)^ (z + Zoo)^, 

as the distortion derivative of an affine structure on C \ {— Zoo, Zoo}. 
It would be nice to: 

• have an isomorphic model of this affine surface as polygonal pieces in 
glued along their edges by affine maps. 
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• precise in which sense it can be considered as a limit of the sequence of 
afline surfaces «S, independently of their conformal uniformization to C \ 

{zi,Z2,Z3,Z4}, 

• reprove the convergence of the corresponding function r], independently 
of the previous sections, i.e. without computing and estimating all these 
integrals; in particular, the uniqueness of the extracted limits of Re (zi) and 
Im(2i)logif would follow from the uniqueness of (normalized) conformal 
uniformization. 

In this section, we will present a candidate isomorphic model for the limit afBne 
surface. 

5.3.2. The candidate. Let us motivate the definition of this candidate. 




Figure 16: Part of the construction of the Riemann surface over C corresponding 
to the developing map of <S with Chart 1 as a germ. The picture illustrates the case 
K = 2. The arrows indicate which sides shall be glued together. 

Given an afhne surface S, the affine charts can be lifted to its universal cover S 
and define an affine surface structure on it. Moreover, since <S is simply connected, 
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any germ of an afRne chart extends uniquely to all of S. This extension is called 
the developing map. It is not necessarily injective anymore. This developing map 
5 — > C can be seen as a structure of "Riemann surface over C" . Do not be mislead 
by the terminology: a Riemann surface over C is much more rigid than a Riemann 
surface; the transition maps between charts are required to be the identity; so it is 
more rigid than an afhne surface or a translation surface. 

For a given K < +00, let us look at a little variant: consider the usual afhne 
surface S, which we recall is defined afhne gluings of two pieces: the complement 
of a square (Chart 1) and a rectangle (Chart 2). Instead of looking at the universal 
cover of S, let us look at the universal cover of 5 U {00}. Remove to this universal 
cover all points that project to 00. Call this afhne surface S' . It has also a well 
defined developing map, because such a map can be dehne by gluing together, by 
z t—f z, afhne images of the square's complement and of the rectangle. This is 
illustrated for if = 2 in hgurefTHl 




Figure 17: A flat model of the limit afhne surface. The arrows indicate which side 
shall be glued together. The spiraling stacks represent infinitely many copies of the 
slit plane C \ (—00, 0], glued together along the slit by z 1-^ z; their afhne structure 
is the one induced by the inclusion in C. 

Playing with this, letting K — > +00, we end up suspecting that the following 
surface is a natural limit of the afhne surface, adapted to the normalization at 00 
(a.k.a. $'(z) = z + + o(l) in Chart 1). It is dehned as follows, and is illustrated 
in figure fTTl Let = C \ Sq. Let Be be defined by "Rez < 1 & |Imz| < 1". It is 
seen as a limit of the rectangle, whose width went to +00 on the left. Similarly, let 
Br be defined by "Rez > — 1 & |Imz| < 1". Glue Be to Oc along the open segment 
{1 — i,l + i) (recall this means that we thicken a little bit the open sets along this 
segment and glue them by z 1-^ z). Glue Br to Oc along (—1 — + i). Glue 
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the upper side and the lower side of the square bounding Oc, ends excluded, by a 
translation. Now we will add four chains of pieces on respectively the upper side 
of S^, its lower side, the upper side of and its lower side. We will describe only 
the first chain, the others are symmetric. Let U2n = {2 G C | Re z > or Im z > O} 
and t/2n+i = {z G C I Rez < and Imz < O}. Glue Uq to Bi along (—00, 0), U2n 
to C/2n+i along (— cx),0) x i and U2n+i to C/2n+2 along (— oo,0). Gluing together 
the Un gives a piece that can also be seen as the subset "arg(z) e (—00, tt)" of the 
universal cover of C* on which we would have put the following affine structure: the 
natural affine structure on C is first pulled back on C* by the inclusion map, then 
pulled-back by the universal cover of C*. This is represented as spirals in figure flTl 
The uniformization to C \ {xoo,— a^oo} is supposed to yield figure [T51 which was 
drawn using the formula for the limit distortion derivative rjao and drawing the 
curves corresponding to the black lines in figure [T71 




Figure 18: The corresponding regions, after uniformization. 
6. Conclusion 

The story is not finished. We may wonder what is the limit when the major axis 
of the ellipses are slanted. The same approach as in this article works, with slight 
complications; the limit is less surprising. We may want to let the value of /ig (= 
the constant value of /i within the square) vary on all of D and even beyond: it 
seems to be possible to do this; one should find an interpretation for it. We may let 
jjL tend to 1 along directions not contained in R and study how the limit depends 
on the direction. 

The fact that when a finite type affine surface can be uniformized onto the 
Riemann sphere § (finitely punctured) the affine coordinates on § are given by a 
generalized Schwarz-Christoffel formula was probably already known, but is not 
so easy to find in the literature^ Concerning the idea of using affine surfaces to 



"'^"it can be found in a slightly different form, like in [G], according to whom it was already 
known from Schwarz and Schlafli. 
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uniforniize Beltrami forms that are piecewise constant on polygons, the author is 
not aware of this having already been tried before. This should be studied further. 

Concerning afiine surfaces and the distortion derivative: if a holomorphic func- 
tion 7] has a singularity of polar type of order > 1, can we find and classify models 
of the afiine surface structure that r] defines near the singularity? Can we do the 
same for projective structures, i.e. with the Schwarzian derivative replacing the 
distortion derivative? 
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